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Summary 

A predictive model for the three-dimensional failure 
behavior of engineering polymers has been developed in a 
recent NASA-sponsored research program. This model 
acknowledges the underlying molecular deformation mech- 
anisms and thus accounts for the effects of different chemical 
compositions, crosslink density, functionality of the curing 
agent, etc., on the complete nonlinear stress-strain response 
including yield. The material parameters required by the model 
can be determined from test-tube quantities of a new resin in 
only a few days. Thus, we can obtain a first-order prediction 
of the applicability of a new resin for an advanced aerospace 
application without synthesizing the large quantities of material 
needed for failure testing. This technology will effect order- 
of-magnitude reductions in the time and expense required to 
develop new engineering polymers. 

Introduction 

The development of models for predicting the important 
engineering properties of aerospace epoxy resins has been the 
focus of a recent research effort at Purdue University 
sponsored by the NASA Lewis Research Center. The objective 
of this research program is to significantly reduce both the 
time and expense required to evaluate new candidate resins 
for high-performance aerospace applications. In this report we 
summarize our research progress to date. 

The development of new polymer matrix resins for use in 
fiber-reinforced composites is both expensive and time con- 
suming. A major expense in the development process is the 
synthesis and subsequent physical characterization of candidate 
resins. Specifically, in most engineering applications the resin’s 
mechanical properties, and in particular its failure properties, 
are of critical importance. In order to even partially evaluate 
the mechanical performance, large (i.e., pound) quantities of 
a new resin are required; however, a polymer chemist typically 
synthesizes gram quantities of a new resin. The scaleup of the 
synthetic scheme to produce large quantities of a new resin 
is usually expensive and time consuming; additionally, there 
is no a priori assurance that the particular candidate resin 
selected for scaleup will have the required physical properties. 
Thus, a repetition of the entire process may be required. 
Moreover, since the deformation and failure properties 
of aerospace resins depend on temperature, deformation rate, 
mode of loading, etc., complete characterization of the 


mechanical properties is quite involved, even when large 
quantities of a new resin can be readily obtained. 

In order to significantly reduce the time and expense required 
to evaluate a new resin for aerospace applications, a new 
constitutive equation has been developed. This constitutive 
equation acknowledges the molecular processes that control 
the deformation and ductile failure in amorphous polymers. 
The only material properties required by the constitutive 
equation are the equilibrium pressure-volume-temperature 
(PVT) response and the linear viscoelastic shear and bulk 
moduli, properties that can be quickly evaluated from test-tube 
quantities of a new polymer resin. As shown schematically 
in figure 1, the nonlinear deformation and failure properties 
can be predicted via the constitutive equation, effecting order- 
of-magnitude reductions in the time and expense required to 
evaluate a new resin. In the next section of this report, we 
discuss briefly the constitutive equation and review some of 
its predictive capabilities. 

The new nonlinear viscoelastic constitutive equation can 
effect a considerable reduction in the time and expense required 
to evaluate a new polymer resin; however, it would be highly 
advantageous if the important engineering properties of the 
resin could be estimated from chemical structure alone. Such 
predictive capabilities would allow chemists who work with 
synthetic polymers to screen a large number of potential 
molecular architectures prior to actual synthesis. Also, a model 
that only requires the chemical structure would provide guidance 
to the polymer chemist, since it is not otherwise obvious what 
types of molecular structures are most effective for enhancing 
important properties such as toughness and fatigue life. Some 
progress has been made in developing group contribution 
methods for estimating the linear viscoelastic properties required 
by the nonlinear constitutive equation. The research results to 
date are summarized in the Group Contribution Methods for 
Linear Viscoelastic Shear Modulus section. 


Symbols 

a shift factor 

af, a\ shift factors for a and /3 mechanisms, respectively 
B material constant 

E Young’s modulus 

E g glassy tensile modulus 

E(t) tensile modulus 
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Figure 1.— Schematic of way in which new nonlinear viscoelastic constitutive equation can significantly reduce time and expense required to evaluate new 
aerospace resins. 


G(t) linear viscoelastic shear modulus 
G'(o3,t) dynamic storage modulus 
G"(o),t) dynamic loss modulus 

H ai , Hpj discrete spectra values for a and 0 mechanisms, 
respectively 

I identity tensor 

K (0 bulk modulus 

P pressure 

P(t) pressure history 

P eq equilibrium pressure at T(t) and P(t) 

S entropy 

S 0 initial entropy 

T absolute temperature 

T g glass transition temperature 

T 0 initial temperature 

7X0 temperature 

T%f,T $ e f reference temperatures for a and 0 mechanisms, 
respectively 

t time 

t* reference time 


V 0 initial volume 

V(t) volume 

e strain 

e(r) displacement 

e strain rate 

£ material time 

a stress 

a(t) stress response 

o) experimental frequency 

Nonlinear Viscoelastic Constitutive 
Equation 

A new constitutive equation for amorphous polymers has 
been developed that can describe (1) the nonlinear viscoelastic 
behavior including yield when the polymer is in the solid state, 
(2) the glass transition phenomenon, and (3) the non- 
Newtonian flow behavior of the polymer melt. Several earlier 
forms of the constitutive equation have been derived intuitively 
(refs. 1 to 3). However, the most general form of the consti- 
tutive equation has been derived (ref. 4 and Lustig, S.R.; Shay, 
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R.M., Jr.; and Caruthers, J.M.: Thermodynamic Constitutive 
Equations With Memory Using a Material Time Scale, to be 
published in Journal of Rheology) via the framework of 
rational mechanics (refs. 5 and 6), where it was assumed that 
the viscoelastic dissipation occurs on a material time scale 
controlled by the nonequilibrium entropy. The material time 
scale incorporates the fact that the rate of viscoelastic relaxation 
changes with the thermodynamic state of the polymer (i.e., 
temperature, pressure, and specific volume). 

Nonlinear Stress-Strain Behavior Including Yield 

Although the complete constitutive equation is quite 
complex, limiting forms are much simpler. For isothermal 
deformations in the solid, the general constitutive equation 
essentially reduces to the following constitutive equation, 
where the stress response a(t) to an arbitrary displacement 
history e(f) is given by reference 2: 
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and G(t) and K (t) are the linear viscoelastic shear and bulk 
.moduli, E(t) is a finite strain tensor computed from the 
displacement strain tensor e(t), t* and £* are material times 
(ref. 6), and a[T(t),V(t)] is a shift function that depends on 
the nonequilibrium entropy S(T,V t P) as suggested by Adams 
and Gibbs (ref. 7). The equilibrium contribution to S was 
determined by using the Simha-Somcynsky equation of state 
(ref. 8). Equation (1) is quite similar to the well-known linear 
viscoelastic constitutive equation, except that the infinitesimal 
strain tensor has been replaced with a finite strain tensor and 
normal time t has been replaced with the material time t * . The 
constitutive equation has no adjustable parameters , and the 
required material properties are (1) the equilibrium PVT 



Figure 2. -Stress-strain curve for diglycidylether of bisphenol-A (DGEBA) 
epoxy resin cured with 4,4-diaminodiphenylsulfone (DDS) and theoretical 
predictions of equation (1). Temperature, 100 °C; stress divided by glassy 
modulus; strain rate, e, 10 -3 sec -1 . 


behavior need to specify the Simha-Somcynsky equation of 
state and (2) the linear viscoelastic shear and bulk moduli. 

Predictions made using equation (1) for the uniaxial exten- 
sion at a constant strain rate of a typical aerospace epoxy resin 
are shown in figure 2, where the PVT and linear viscoelastic 
material properties were measured in independent experiments. 
By examining figure 2, we see that the constitutive equation 
predicted the nonlinear stress-strain curve including yield and 
that the predictions agreed quantitatively with the experi- 
mentally measured stress-strain response. The discrepancy 
between theory and experiment after yield was to be expected 
because the predictions assumed a spatially homogeneous 
deformation whereas the experimental specimen exhibited 
necking. In order to predict spatially nonhomogeneous defor- 
mations such as necking, equation (1) must be combined with 
the equation of continuity and conservation of momentum in a 
finite-element program. Since the yield data were not used to 
set any parameter in the constitutive equation, the quantitative 
agreement shown in figure 2 is a definitive verification of the 
new constitutive equation. In addition to quantitatively 
predicting the yield phenomenon, the constitutive equation is 
able to predict the effects of strain rate, temperature, and 
multiaxial deformations on the yield process (refs. 2 and 3). 

Glass Transition Phenomenon 

Another limiting form of the general constitutive equation 
(ref. 4 and Lustig, S.R.; Shay, R.M., Jr.; and Caruthers, J.M.: 
Thermodynamic Constitutive Equations With Memory Using 
a Material Time Scale, to be published in Journal of Rheology) 
is the PVT behavior in the glass transition region. The time- 
dependent volume response V(t) to an arbitrary temperature 
T(t) and pressure P(t) history is given by 
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P(t)=P eq [T(t),V{t)] + ^^\+h u jj 

( 4 ) 

Jo 3 « 

where e&. = [K(r) /K 0 ] 2/3 - 1 is the voluminal strain, V 0 is the 
initial volume, and P eq is the equilibrium pressure at T(t) and 
V(t) as given by the equation of state. Here again, there are 
no adjustable parameters and the only material properties are 
the equilibrium PVT behavior and the linear viscoelastic bulk 
modulus. The predictions of equation (4) for poly (vinyl acetate), 
as it was cooled from an equilibrium rubber, are shown in 
figure 3. The constitutive equation predicted a discontinuity in 
the slope of the volume-temperature curve that would be 
identified with the glass transition temperature T g . For 
poly (vinyl acetate) T g is 32 °C, in quantitative agreement with 
the theoretical prediction. Since T g was not used to set any 
model parameter, the quantitative agreement shown in figure 3 
is a definitive test of the constitutive equation’s ability to describe 
the glass transition phenomenon. Equation (4) also accurately 
predicts the effects of the cooling rate on the glass transition. 

Fluid Behavior 

The general constitutive equation can be simplified for the 
case of steady shear of an uncrosslinked polymer fluid, although 
the form is not as simple as equations (1) and (4) and thus 
will not be shown here. This constitutive equation again has 
no adjustable parameters and only requires the equilibrium 
PVT response and the linear viscoelastic shear and bulk 
moduli. Preliminary prediction of the non-Newtonian shear 
viscosity for a model poly(vinyl acetate) material is shown in 
figure 4. The constitutive equation predicted a zero shear 
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Figure 3.— Specific volume versus temperature predicted by equation (4) for 
poly(vinyl acetate) cooled at 5 deg C/hr. Pressure, 1 bar. 
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LOG SHEAR RATE 

Figure 4.— Predictions of nonlinear viscoelastic constitutive equation of strain- 
rate-dependent shear viscosity for single-relaxation-time idealized poly(vinyl 
acetate) material at 80 °C. 

viscosity at low shear rates and shear thinning viscosity at 
higher shear rates. The shear thinning region could be described 
by a power law with a slope of approximately -1/2, which 
was in qualitative agreement with experimental data. The 
constitutive equation also predicted a primary normal stress 
that was also shear thinning and in qualitative agreement with 
experimental data. Although much work remains before the 
non-Newtonian fluid predictions of the constitutive equation 
are fully evaluated, the initial results shown in figure 4 look 
quite promising. 

Group Contribution Methods for Linear 
Viscoelastic Shear Modulus 

An extensive experimental study of the linear viscoelastic 
shear modulus was performed for a series of diglycidylether 
of bisphenol-A (DGEBA) epoxy resins cured with 
4,4-diaminodiphenylsulfone (DDS). The dynamic linear 
viscoelastic storage and loss moduli, G' and G" , were 
measured as functions of frequency from liquid nitrogen 
temperatures to well above T g (ref. 9). The complete linear 
viscoelastic response is usually determined by time-temperature 
superposition of the frequency- or time-dependent isotherms 
(ref 10); however, time-temperature superposition is not 
applicable for the DGEBA/DDS materials, since they exhibit 
multiple relaxation mechanisms. 

A novel procedure has been developed (ref. 1 1) to extract 
the linear viscoelastic shear response when there are two 
relaxation mechanisms: the a mechanism in the glass-to-rubber 
transition and a sub-7^ j8 mechanism. We assumed that the 
viscoelastic response of the material was a linear combination 
of the a and (3 response; that is, 

G>,r) = G' a (a^,7“ f ) + Gp (a^Ji f ) (5) 

and 


G'(«.D = c; (afw.rSf) + g ; (af«,7* f ) (6) 

where and 7"" ef ,rf ef are the shift factors and reference 
temperatures for the a and 0 mechanisms, respectively; T is 
the experimental temperature; and co is the experimental 
frequency. It is implicitly assumed in equations (5) and (6) 
that the relaxation spectra of the a and 0 mechanisms are 
additive but that the time-temperature shifts are different. For 
discrete spectra of a and 0 relaxation times, the G"(cu) 
isotherms are given by 


G"(u,T)=Yj GZfa.T) + £ GZj(u,T) 

all / all j 


all i 


0)0 T T i , Y* M ua^Tj 

1 + (c*iS^) 2 L ^ 1 +(wflfry) 2 


(7) 


where // af and are the discrete spectra values for the a 
and 0 mechanisms. The G"(o),T) response is experimentally 
measured at a variety of frequencies and temperatures. 
Subsequently, equation (7) is written in a matrix form, where 
H ai and H@j are elements of an unknown column vector H . 
By using the af and a $ shift functions determined from 


shifting the G" data, the vector H can be determined from 
the experimental G" isotherms. 

Because the set of equations is nearly singular, the conven- 
tional least-squares solution procedure produces oscillatory 
spectra. The technique of regularization with quadratic pro- 
gramming (RQP) is an appropriate method for solving ill-posed 
problems such as equation (7) (ref. 10). We have used the RQP 
method with proper scaling of the variables in equation (3) 
to accurately determine H ai and H^. Details of the method 
are presented elsewhere (ref. 9). The log a“and log af shift 
factors are shown in figure 5 for the series of epoxy resins. 
The individual G" responses, as determined from the spectra 
for the a and 0 mechanisms, are shown in figure 6 for the 
series of epoxy resins. 

An examination of the log a j data in figure 5 shows that the 
a-mechanism shift factor is the same for all the epoxy materials 
when referenced to T g . The difference in log af below T g is 
a consequence of different thermal histories during formation 
of the glass rather than an effect of the chemical structure. The 
G" response for the a-mechanism is essentially identical for 
all the epoxy resins, with the exception of the two most tightly 
crosslinked materials, which have greater strength at high 
frequencies. Further, the G" response for the 0 mechanism 
was essentially the same for all the epoxy resins investigated. 
The data shown in figures 5 and 6 are the beginnings of a group 
contribution method for the linear viscoelastic shear modulus 
of DGEBA epoxy resins. Group methods for other diepoxide 
systems have also been investigated (ref. 12). 




(a) a mechanism. Reference temperature, T g . 

(b) 0 mechanism. Reference temperature, -70 °C. 

Figure 5.— Independent log a“ and log of shift factors for series of DGEBA/DDS epoxies. 
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EXPERIMENTAL FREQUENCY, log g, rad/sec 

(a) a mechanism. 

(b) (3 mechanism. 

Figure 6.— Loss modulus response for a and (3 mechanisms for series of DGEBA/DDS epoxies. 


Summary of Results 

From this study to develop models for predicting the 
important engineering properties of aerospace epoxy resins, 
the following results have been noted: 

1 . A new constitutive equation has been developed that can 
describe (a) the nonlinear viscoelastic behavior including yield 
when the polymer is in the solid state, (b) the glass transition 
phenomenon, and (c) the non-Newtonian flow behavior of the 
polymer melt. 

2. The new constitutive equation has no adjustable 
parameters and the required material properties are (a) the 
equilibrium pressure-volume-temperature behavior needed to 
specify the Simha-Somcynsky equation of state, and (b) the 
linear viscoelastic shear and bulk moduli. 

3. The constitutive equation predictions for a typical epoxy 
resin agreed quantitatively with the experimentally measured 
stress-strain response. 

4. The specific volume versus temperature response for 
poly(vinyl acetate) agreed quantitatively with measured data 
and accurately predicted the glass transition temperatures T g . 

5. The constitutive equation accurately predicted the strain- 
rate-dependent shear viscosity for single-relaxation-time 
poly (vinyl acetate) material. 

6. A novel procedure has been developed to extract the linear 
viscoelastic shear response when there are two relaxation 
mechanisms (a and sub -T g ($) operating. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, March 30, 1989 
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